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Abstract. We discuss different parameterizations of the C^{q 2 ) NA axial form factor, fitted to the old Argonne bubble 
chamber data for pion production by neutrinos, and we use coherent pion production to test their low q 2 behavior. We find 
moderate effects that will be difficult to observe with the accuracy of present experiments. We also discuss the use of the Rein- 
Sehgal model for low energy coherent pion production. By comparison to a microscopic calculation, we show the weaknesses 
of some of the approximations in that model that lead to very large cross sections as well as to the wrong shapes for differential 
ones. Finally we show that models based on the partial conservation of the axial current hypothesis are not fully reliable for 
differential cross sections that depend on the angle formed by the pion and the incident neutrino. 
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INTRODUCTION 

A proper understanding of coherent pion production is 
very important in the analysis of neutrino oscillation ex- 
periments since pion production is a source of back- 
ground [ 1 ] . With pions being mainly produced through 
the excitation of nucleon resonances, coherent produc- 
tion can be used to extract information on axial form fac- 
tors for the nucleon-to-resonance transition. 

In coherent production the nucleus remains in its 
ground state and the reaction is controlled by the nucleus 
form factor. The nucleus form factor favors small val- 
ues of the nucleus momentum transfer squared t. Small 
t values imply in this reaction small q 2 (square of the 
lepton momentum transfer). For small q 2 , coherent pion 
production is dominated by the divergence of the axial 
current and it can thus be related to the pion-nucleus co- 
herent scattering through the partial conservation of the 
axial current (PCAC) hypothesis. 

Experimental analyses of the coherent reaction rely 
on the Rein-Sehgal (RS) model [0] which is based on 
PCAC. In the RS model the pion-nucleus coherent cross 
section is written in terms of the pion-nucleon elastic 
cross section by means of approximations that are valid 
for high neutrino energies and small t and q 2 values. As 
pointed out in Refs. [3, 4], those approximations are not 
reliable for neutrino energies below/around 1 GeV and 
light nuclei, like carbon or oxygen. 

There are other approaches to coherent production that 
do not rely on PCAC but on microscopic models for 
pion production at the nucleon level yl @]. The domi- 



nant contribution to the elementary amplitude at low en- 
ergies is given by the A-pole mechanism (A excitation 
and its subsequent decay into pion nucleon). Medium 
effects like the modification of the A mass and width 
in the medium, final pion distortion, evaluated by solv- 
ing the Klein-Gordon equation for a pion-nucleus op- 
tical potential, as well as nonlocalities in the pion mo- 
mentum, are very important and are taken into account 
in microscopic calculations. In the microscopic model of 
Ref. D6D background terms were included on top of the 
dominant A-pole contribution. The least known ingredi- 
ents of the model are the axial NA transition form factors, 
of which C5 gives the largest contribution (See Eq. 1 of 
Ref. ||7|] for a form factor decomposition of the A^A weak 
current). Besides, within the Adler model [8] used in 
Ref. 1 6], Cj determines all other axial form factors. This 
strongly suggested the readjustment of that form factor 
to the experimental data, which the authors did by fit ting 
the flux-averaged v^p — > pL~ pK + Argonne (ANL) §,[11 
^-differential cross section for pion-nucleon invariant 
masses W < 1.4 GeV, for which the model should be 
appropriate. They found C5 (0) = 0.867 which is some 
30% smaller than the value predicted by the off-diagonal 
Goldberger-Treiman (GT) relation that predicts Cj (0) = 
1.2. Background terms turn out to play a minor role in 
coherent production [3] where a reduced C5 (0) value 
gives rise to smaller cross sections. Coherent production 
is dominated by the axial current which in microscopic 
models is, in its turn, dominated by the C5 axial form 
factor. Coherent production could then be used to test the 
validity of different C^(q 2 ) dependences in the small q 2 



region. 

In this contribution we will concentrate on the two is- 
sues mentioned before: first, we will try to see how sen- 
sitive coherent production is to different C^(q 2 ) param- 
eterizations proposed in the literature. Second, we will 
present results that show how and why the RS model fails 
for low energy pion coherent production on light nuclei. 



C^{q 2 ) IN COHERENT PRODUCTION 

As mentioned before, in Ref. Jfl the authors made a fit of 
C5 (q 2 ) to ANL data. For C5 they took the q 2 dependence 
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of Ref. HI 
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and from the fit they obtained 

Cs(0) = 0.867 ± 0.075, M AA = 0.985 ±0.082 GeV. (2) 

This fitted axial mass in the weak NA vertex is in good 
agreement with the estimates of about 0.95 GeV and 0.84 
GeV given in the original ANL reference II 1 0TI and in the 
work of Ref. Il2ll . On the other hand, a correction of the 
order of 30% to the off-diagonal GT relation value was 
found for CH0). As shown in Ref. Q6D the agreement with 
ANL total cross sections improved with the fitted values 
for C5 (0) and M A a- On the other hand it is also shown 
that Brookhaven (BNL) bubble chamber data II 1311 and 
ANL data are not fully compatible and that BNL data 
alone would favor a C5 (0) ~ 1 .2 value as given by the 
off-diagonal GT relation. 

A different approach has been followed by Leitner 
et al. in Ref. lfl4fl . There, the authors use a different 
parameterization for (q 2 ) 
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aq 
b — q 



■)/{l-q 2 /M AA 2 



(3) 



in which C5 (0) is kept to its off-diagonal GT relation 
value C5 (0) = 1 .2, while a, b and M A a are treated as free 
parameters. One can accommodate a larger C5 (0) value 
by changing the q 2 dependence. In fact very small — q 2 
values are not very relevant due to phase space and what 

value in the region 

-<t 

the A-pole term alone they got Hf 

a = -0.25, b = 0.04 GeV 2 , M AA = 0.95 GeV. (4) 

But background terms are important at the nucleon level 
and they should be included in the calculation. A new fit 



is actually important is the C5 ( 

around -a 2 ~ 0. 1 GeV 2 . When fitting the ANL data with 
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FIGURE 1. Different C^(q l ) parameterizations fitted to 
ANL data. 



including background terms leads to 

a = -0.3861 ±0.198, b = 0.01536 ±0.0310GeV 2 , 
Maa — 0.952 ± 0.205 GeV. (5) 

Both the fit in Eqs. Q] and [2] and the one in Eqs. [3] and 
[5] give a good description of ANL data even though 
they use quite different values for Cj (0). The different 
q 2 dependence compensates for the initial difference at 
q 2 = 0. Note, however, the large errors in Eq.[5]that point 
at large correlations between the different parameters. In 
Fig.Q]we compare the two form factor parameterizations. 
The main differences are in the -q 2 < 0.025 GeV region. 
The corresponding axial radii are R\ =0.56+ 

(109 f m f° r 

the parameterization given by Eqs. Q] and |2j and R\ = 
6.4^g ^fm 2 for the one given by Eqs.[3]and[5] The large 
negative error in the latter case is a reflection of the large 
statistical errors in Eq. [5] Both things point to the fact 
that Eq.[3]might not be a good parameterization. 

In the following we will try to see if coherent pion 
production can give us extra information on the validity 
of the two different C^(q 2 ) discussed above. Here we 
shall use the coherent production model of Ref. 13J . 

In Fig. [2] we show the q 2 differential cross section for 
charged current (CC) and neutral current (NC) coherent 
production on carbon and for an incident neutrino energy 
E = 600 Me V. In both cases one can see differences for 
low — q 2 due to the different C^(q 2 ) form factor used. 
These effects are more relevant in the NC case where 
smaller — q 2 values are not suppressed by phase space. 



As 



da 
dq- 



can not be measured for the NC channel one has 



to rely on total cross sections. There, the change is a 
mere 3.4% for the CC channel, while for the NC process, 
where lower — q 2 can be reached, the change is 16.6%. 

In Fig. [3] we show the NC differential cross section 
with respect to cos Q K , being Q K the angle formed by the 
pion and the incident neutrino. In this case there is no 
difference in shape as can be appreciated by comparing 
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FIGURE 2. 



da 



differential cross sections for CC and NC 



neutrino induced one-pion production on carbon obtained with 
different C^(q 2 ) parameterizations. The incident neutrino en- 
ergy is E = 600 Me V. 



the solid line and the up triangle line. A similar result is 
obtained for the CC case and for the NC/CC differential 
cross section with respect to the pion kinetic energy. 

Finally in Fig. [4] we show total cross sections. For the 
CC case a cut \k^ \ > 450 Me V on the final muon momen- 
tum, as the one imposed by the K2K Collaboration 01511 . 
has been applied. The changes are more significant for 
the NC case. In Table [T] we show the CC results con- 
volved with the K2K flux and the NC ones convolved 
with the MiniBooNE flux. Details of the convolution can 
be found in Ref. [3]. We also show for comparison the 
experimental data. The change in the parameterization 
for C$(q 2 ) amounts to an small increase of 3% for the 
CC/K2K case. Both results are below the present experi- 
mental upper bound. For the NC/MiniBooNE case there 
is an increase of 10.7%, with the results of both parame- 
terizations being a factor of 2 smaller than present exper- 
imental data IU6I1 . 

Coherent pion productions shows a moderate depen- 
dence on different C^(q 2 ) parameterizations fitted to the 
ANL data. With present experimental accuracy it can not 
be used to constraint the form factor small q 2 behavior. 
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FIGURE 3. -rS&ir differential cross section for the v 12 C -> 

Li COS 

v 12 C 7t a reaction obtained with different C^(q 2 ) parameteri- 
zations. The incident neutrino energy is E = 600 MeV 
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FIGURE 4. Total cross sections, as a function of the incident 
neutrino energy, for CC and NC neutrino induced one-pion 
production on carbon obtained with different C^{q 2 ) param- 
eterizations. For the CC case a cut [k^ | > 450 MeV on the final 
muon momentum has been applied. 



Faxis, and calling (j) Kq to the pion azimuthal angle in the 
XYZ frame, Lorentz invariance allows us to write 



REIN-SEHGAL MODEL FOR LOW 
ENERGY COHERENT PRODUCTION 

The RS model is based on PCAC, so it is worth moti- 
vating how one constructs a general PCAC based model. 
Let us look, for simplicity, at the NC process 

v,(k(E,k)) + ,yV gs -> v, (k'(E',k')) + jK gs + 7Z°(k K ) . (6) 

Defining q = k — k', y = q°/E, t = (q — k % ) 2 , taking q 
along the positive Z axis and k x k along the positive 
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where G is the Fermi decay constant, K = q° + j-^, 
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Besides, 



with the nucleus mass, and u.v — — ^et 
<7r.l,s stand for cross sections for right, left and scalar 
polarized intermediate vector mesons, si is not a proper 
cross section and it contains all the dependence on ty nq . 
As shown in Ref. [ 4] , Eq. |7]should be the starting point to 



TABLE 1. Total cross sections for CC and NC neutrino induced one-pion production on carbon 
convolved with the K2K and MiniBooNE fluxes. Different C^(q 2 ) parameterizations have been used. 
Experimental results are also shown. 
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CC V M + 12 C 
CC v M + 12 C 


K2K 
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Eqs.[3and[2] 
Eqs.[3]and[5] 


4.68 
4.82 


J < 7.7 [15] 


NC v M + 12 C 
NC v M + 12 C 


MiniBooNE 
MiniBooNE 


Eqs.[Hand|2] 
Eqs.[3]and[5] 


1.97 
2.18 


| 4.54±0.04±0.71 [16] 



evaluate cross sections with respect to 9 n . PC AC based 
models take instead as a starting point 
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which is obtained from Eq. [7] after integration on ty aq , 
and they further assume 
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The latter is incorrect for q 2 ^ and it will have con- 
sequences when evaluating cross sections with respect to 
variables that depend on Q K . 

For q 2 = only (Js contributes and one finds that 

is given as the modulus square of the hadronic 
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matrix element of the divergence of the weak current. 
Since the vector NC current is conserved one is left with 
the divergence of the axial current that can be related, 
through PCAC, to the pion-nucleus elastic cross section 
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with f n = 92. 4 Me V. Neglecting the nucleus recoil (k® — 
q°) and including a form factor Ga = 1/(1 — q 2 /m\) for 



q ^ , one arrives at 
da 
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This is the Berger-Sehgal model for 7T° coherent produc- 
tion lfl7ll . In the RS model they further approximate 
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with 7vi(f ) the nucleus form factor, F a b s a t— independent 
eikonal absorption factor that takes into account the dis- 
tortion of the final pion, and da{ - n N ~> the elastic pion- 
nucleon differential cross section at t = 0. 



It is worth modifying the RS model by improving 
some of its approximations |4]. The t = approximation 
can be eliminated altogether by substituting 
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where nsp stands for the non-spin-flip part of the pion- 
nucleon elastic cross section. Besides, the pion distortion 
can be improved, still in an eikonal approach, with the 
replacement 
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Finally one can use ^=fp instead of ^-y-. 

In Fig.[5]we evaluate J^j for NC coherent production 
on carbon at E = 0.5 GeV. In order to check the t = ap- 
proximation of the RS model, no distortion is included in 
the calculation. We compare with the modified RS model 
where the t = approximation has been removed. We 
also compare with the microscopic calculation of Ref. Ot] 
that we consider a good model for coherent production at 
these energies. For simplicity, in the microscopic calcula- 
tion we have only kept the dominant A contribution and, 
in order to make the comparison meaningful, we have 
taken Cj (0) = 1.2 (to fix normalization) and we have not 
included any in medium correction for the A or any pion 
distortion. We see that the modified RS model compares 
very well with the microscopic calculation, whereas the 
original RS model fails to reproduce both the size and 
the shape of the differential cross section. 

In the same conditions as before, we show in Fig. [6] 
the ^ differential cross section on carbon at E = 1 GeV. 
Once more, we see the size and shape of the RS calcula- 
tion is very different from the microscopic one, showing 
the inadequacy of the t = approximation. In this case 
the modified RS is also in disagreement with the micro- 
scopic calculation. This is mainly due to the approxima- 
tion in any PCAC based model encoded in Eq.|9]and, to a 
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FIGURE 5. 4^r differential cross sections for the v 12 C -» 
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V 12 C 7T° reaction obtained within different models. The in- 
cident neutrino energy is E = 0.5 GeV. No distortion or in 
medium corrections were considered in the calculations. 
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v 12 C 71° reaction obtained within different models. The inci- 
dent neutrino energy is E = 1 GeV. No distortion or in medium 
corrections were considered in the calculations. 



lesser extent, to the neglect by PCAC models of the (7r,l 
contributions. This assertion can be checked by eliminat- 
ing from the microscopic calculation the non-PCAC con- 
tributions, i.e. taking only (7$ and assuming Eq. [9] The 
results are shown in Fig.|6]with up triangles. We now see 
the agreement with the modified RS calculation is good. 
This shows in fact that PCAC based models are not fully 
reliable to give differential cross sections with respect to 
variables that depend on Q K . 

Last, in Fig. [7] we present results for ^ but now with 
distortion and in medium corrections. Again, the sizes 
and shapes of the RS model predictions are very differ- 
ent from the microscopic ones. The same is true for the 
modified RS calculation, even when compared to the mi- 
croscopic one with only PCAC contributions. While the 
shapes of the differential cross sections of the two lat- 
ter calculations are similar, the integrated cross sections 
differ by a factor « 1.7. This hints at the inadequacy of 
eikonal factors to account for the final pion distortion or 
in medium modifications of the 7CN cross section. We see 
again that the approximations in PCAC models give rise 
to the wrong distributions for variables depending on n . 

As shown, the t = approximation and the eikonal 
distortion used in the RS model are not adequate for low 
energy coherent production on light nuclei. Besides, and 
due to the assumption in Eq. [9] any PCAC model would 
predict the wrong shape for differential cross sections 
with respect to variables that depend on 6 n . 
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